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The HIJING (Heavy-ion Jet Interaction Generator) Monte Carlo model is updated with the latest 
parton distributions functions (PDF) and new set of the parameters in the two-component mini- 
jet model that controls total p + p cross section and the central pseudorapity density. We study 
hadron spectra and multiplicity distributions using the HIJING 2.0 model and compare to recent 
experimental data from p + p collisions at the LHC energies. We also give predictions of hadron 
production in p + p, p + Pb and Pb + Pb collisions at the full LHC energy. 
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I. INTRODUCTION 

The HIJING (Heavy Ion Jet INteraction Generator) 
Monte Carlo model was developed to study hadron pro- 
duction in high-energy nucleon-nucleon, nucleon-nucleus 
and nucleus- nucleus collisions [lHH- It was based on a 
two-component geometrical model of mini-jet production 
and soft interaction and has incorporated nuclear effects 
such as nuclear modification of the parton distribution 
functions and jet quenching via final state jet medium 
interaction. It can reproduce most features of hadron 
production in p + p{p) and p + A collisions up to the Fer- 
milab Tevatron energies. With some modification of the 
string configuration in the soft sector of particle produc- 
tion, it can also reproduce the bulk hadron spectra and 
the approximately the suppression of high px hadrons 
due to jet quenching in central rapidity region of A + A 
collisions up to the RHIC energies [J-Q- ^ has been 
widely used to simulate hadron production in p + A and 
A + A collisions for designs of new detector systems and 
provide initial conditions forparton and hadron cascade 
models such AMPT model 0. 

The core component of HIJING is the two-component 
model for beam parton interaction which was proposed 
to model the energy dependence of total cross section^- 
[T3 | and particle production p^H18| in high energy hadron 
collisions. The two-component model was also extended 
fl9l . [20( | and later incorporated in the HIJING model to 
describe initial parton production in high-energy heavy- 
ion collisions. In this two-component model, one assumes 
that parton interaction in high-energy nucleon-nucleon 
collisions can be divided into soft interaction and hard 
or semi-hard interaction with jet production. A cut-off 
scale po in the transverse momentum of the final jet pro- 
duction has to be introduced below which the interaction 
is considered non-perturbative and can only be charac- 
terized by a finite soft parton cross section a so ft- For 
jet production with transverse momentum px > Po, the 
cross section and jet spectrum are assumed to be given 
by perturbativc QCD parton model. Jet cross sections in 



collinear factorized perturbative QCD in turn depend on 
the parton distribution functions (PDF) that are param- 
eterized from a global fit to the available experimental 
data of deeply inelastic scattering (DIS) of lepton and 
nucleon, Drell-Yan lepton pair, direct photon and jet pro- 
duction in p + p{p) collisions. 

The two parameters <J S oft and po in HIJING are de- 
termined phcnomcnologically by fitting the experimental 
data of total cross sections and hadron multiplicity in 
p + p/p collisions, which should depend on the param- 
eterization of nucleon PDF's. In the original version of 
HIJING, Duke-Own parameterization [21| of PDF's arc 
used to calculate the jet production cross section with 
Pt > Po which are adequate for description of jet pro- 
duction at RHIC and maybe Tevatron energies. How- 
ever, for p + p collisions at much higher energies such as 
the LHC energies, jet production processes involve initial 
beam partons with fraction momentum x ~ 10~ 4 where 
Duke-Owen parameterization of PDF's is no longer valid. 
More modern paramcterizations should be used. In this 
progress report, we introduce the HIJING 2.0 model in 
which we use more modern set of parameterized PDF's. 
We have to perform a global fit of the total cross sections 
and hadron multiplicities in p + p(p) collisions to deter- 
mine the two parameters, a S ojt and po, and their energy 
dependence. We will use the new HIJING 2.0 model to 
study hadron production in p+p at the LHC energies and 
compare to the recently published experimental data. We 
also provide predictions of hadron multiplicities in p + p 
and Pb + Pb collisions at the full LHC energies. 



II. TWO-COMPONENT MODEL IN HIJING 

In the two-component model [D-d, |8H1~I] incorporated 
in HIJING, one assumes that events of nucleon-nucleon 
collisions at high energy can be divided into soft and 
hard processes with at least on pair of jet production 
with px > Po- The inclusive jet cross section Oj et in the 
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depends on the parton-parton cross section a ab and par- 
ton distribution functions fa{x,Pr) in a parameterized 
form of experimental data of other high-energy DIS and 
nucleon-nucleon collisions, where the summation runs 
over all parton species, y\ and 7/2 are the rapidities of the 
scattered partons, x\ and x 2 are the light-cone momen- 
tum fractions carried by the initial partons, and K ss 2 
accounts for the next-to-leading order (NLO) corrections 
to the leading order jet cross section. 

Within the eikonal formalism [8l-[i~3l. [23j . the nucleon- 
nucleon cross sections can be expressed in the impact- 
parameter representation as 
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in the limit that the real part of the parton-parton scat- 
tering amplitude can be neglected. The eikonal func- 
tions s) at an impact parameter b is therefore real 
and can be expressed in term of inclusive jet cross sec- 
tion in pQCD and an effective cross section for the non- 
perturbative soft parton-parton collisions within the two- 
component model [l[, 

X(M) = Xs{b,s) + Xh(b,s) 

= 1 -[a soft T NN (b)+a jet T NN (b)], (6) 

where Tjv^r(6) is the nucleon-nucleon overlap function 
and will be assumed to take the form of the Fourier trans- 
formation of a dipole form factor in HIJING. 

Within the above eikonal implementation of the two- 
component model, one can also calculate the nucleon- 
nucleon cross section [JJJ for no jet and j > 1 number of 
jet production with pt > po, 
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Their sum gives rise to the total inelastic cross section 
in Eq.@. The above eikonal formalism is the base for 
the Monte Carlo simulation of multiple jet production in 
p + p, p + A and A + A collisions. Both cr S oft an d the 



transverse momentum cut-off po are considered param- 
eters in the HIJING model and are fit to the total and 
inelastic nucleon-nucleon cross sections and the hadron 
multiplicity density in the middle rapidity y = at each 
colliding energy. 

One can find the detailed description of HIJING model 
for hadron production in p+p, p + A and A + A collisions 
in Refs. [TJ. The main features of HIJING are: 

(1) Multiple mini-jet production arc simulated accord- 
ing to the above eikonal formalism for each nucleon- 
nucleon collisions at given impact parameter b. The 
kinematics of each pair of jets and the associated initial 
and final state radiation arc simulated using PYTHIA 
modelflTjj. 

(2) Events without jet production (with pr > po) and 
the underlying soft parton interaction in events with jet 
production are modeled by excitation of quark-diquark 
strings with gluon kinks along the lines of the FRITIOF 
model and the DPM modelOU!]. In addition, mul- 
tiple low-pr exchanges among the end point constituents 
arc included. 

(4) A set of impact-parameter-dependent parton distri- 
bution functions is used to include nuclear modification 
of the parton distribution functions inside nuclei. 

(5) A simple model for jet quenching is used to study 
the effect of jet medium interaction in A + A collisions 

M. 



III. UPDATES IN HIJING 

In the default setting of HIJING 1.0 [TJ, Duke-Owens 
parameterization 2l| of PDFs in nucleons is used. With 
Duke-Owens parameterization of PDFs, an energy inde- 
pendent cut-off scale po =2 GeV/c and a constant soft 
parton cross section a so f t = 57 mb are sufficient to repro- 
duce the experimental data on total and inelastic cross 
sections and the hadron central rapidity density in p+p/p 
collisions [2?J ■ Duke-Owens parameterization of PDFs is 
known to be very outdated and one needs to use more 
modern parameterizations from new global fit to experi- 
mental data, especially at the LHC energies when mini- 
jet production reaches to very small- a; region of the par- 
ton distribution where gluon distribution is much higher 
than Dukc-Owcns parameterization. 

Furthermore, with a constant transverse momentum 
cut-off po = 2 GeV/c in HIJING 1.0, the total number of 
min-jets per unit transverse area could exceed the limit 



T AA (b)<T jet p| 
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for independent multiple jet production even in central 
p+p collisions for sufficiently large inclusive jet cross sec- 
tion at high colliding energies, where T A A{b) is the over- 
lap function of A + A collisions and ir/p 2 , is the intrinsic 
transverse size of a mini-jet with transverse momentum 
Po- Therefore, one inevitably has to increase the value of 
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the transverse momentum cut-off po to ensure the appli- 
cability of the underlying two-component model of inde- 
pendent multiple jet production in the HIJING model. 

In the updated version of HIJING 2.0, we will use the 
Gluck-Reya-Vogt(GRV) parameterization 28] of PDFs, 
among many modern parameterizations of the PDFs that 
are available. The gluon distributions in this new param- 
eterization is much higher than the old Duke-Owens pa- 
rameterization at small x and therefore give much larger 
inclusive jet cross section at high colliding energies with 
a fixed value of cut-off p^. One therefore can no longer 
fit the experimental p + p/p data on total and inelastic 
cross sections using a constant cut-off po and the soft par- 
ton cross section a so ft within the two-component model. 
One has to assume an energy-dependent cut-off po(y^s) 
and soft cross section <J S oft(\fs) [29j . Fitting the experi- 
mental values of the total and inelastic cross sections of 
p+p/(p) collisions including those extracted from cosmic 
experimental and the hadron central rapidity density, we 
have the following parameterized energy-dependence of 
the cut-off and the soft parton cross section used in the 
two-component model of HIJING 2.0: 

po = 2.62 - 1.0841og( v / i) +0.299k>g i (y/s) 

-0.02921og 3 ( v / i) + 0.001511og 4 (^), (10) 
a soft = 55.316 - 4.11261og(Vs) + 0.8541og 2 (V"S) 

-0.03071og 3 ( v / s) + 0.003281og 4 (^), (11) 

where the colliding energy ^/s in center-of-mass frame is 
in units of GeV. Shown in Fig. [T] are the calculated to- 
tal and inelastic cross sections using both HIJING 1.0 
and HIJING 2.0 as compared to the experimental data. 
The total inclusive jet cross section and non-pcrturbative 
soft parton cross sections are also plotted for illustration. 
With a constant cut-off po — 2 GcV/c and soft parton 
cross section cr so f t at high colliding energies, HIJING 1.0 
already gives larger total cross section than the cosmic 
data indicate even with the Duke-Owens parameteriza- 
tion of PDFs. With much higher gluon distribution at 
small x in the GRV parameterization used in HIJING 
2.0, one has to introduce a cut-off po and the soft par- 
ton cross section a so ft that increase with colliding energy 
in order to fit the experimental data on the total cross 
section. There are, however, some freedom in fixing the 
values of po and a so ft , which is further constrained by the 
energy-dependence of the central rapidity density of the 
charged multiplicities, as shown in Fig. [2] The increas- 
ing cut-off as required by the experimental data indicates 
that multiple mini-jet production below such cut-off are 
no longer independent and coherent interaction becomes 
important. This might be taken as an indirect evidence 
of gluon saturation at very small x inside a proton in 
proton-proton collisions at very high energies, especially 
at the LHC energies. An alternative approach to effec- 
tively take into account of such gluon saturation is to 
increase the string tension of soft hadron production as 
proposed in Ref. [30|. We choose to focus on the change 
of minijet production in HIJING2.0. 
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FIG. 1: (color online) Total, soft and jet production cross 
sections of pp and pp collisions. The histogram in left panel 
are calculated using HIJING 1.0, while the lower using HI- 
JING2.0. The data are from [Umll . 



We also show in Fig.[3]the transverse momentum spec- 
tra calculated with HIJING 2.0 at different colliding en- 
ergies as compared to the experimental data. HIJING 2.0 
results are all in good agreement with the experimental 
data. 

Note that both HIJING 2.0 calculations and data 
shown in Fig. [2] are for non-single-diffractive (NSD) 
events. The definitions of NSD trigger arc different in 
different experiments leading to different central pseudo- 
rapidity densities in these NSD events. The NSD triggers 
in UA5 experiment require at least one charged parti- 
cle simultaneously in each of the pseudorapidity regions 
at both ends covering 2 < \r/\ < 5.6, while Collider- 
Detector at Fermilab (CDF) NSD events are triggered 
in 3.2 < \r)\ < 5.9. The increase of the central pseu- 
dorapidity density with energy can be attributed to the 
increased mini-jet production in high colliding energies. 
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FIG. 2: (color online) Pseudorapidity distributions of charged 
particles in non-singie-diffractive pp at y/s = 53 GeV, pp col- 
lisions at y/s = 200, 540, 900, 1800 GeV as compared to 
experimental data [33, [33] • 
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FIG. 3: (color online) Invariant inclusive cross sections of 
charged particles in pp at y/s = 53 GeV, pp collisions at y/s = 
200, 500, 900, 1800 GeV. The data come from ^3^4^. Both 
the calculation and experimental data are obtained in the 
region of \rj\ < 2.5 for y/s = 200, 540, 900 GeV, and |r/| < 1.0 
for y/s = 53 and 1800 GeV. 



IV. HADRON SPECTRA IN p + p AND A + A 
COLLISIONS AT THE LHC ENERGIES 

With the updated HIJING 2.0, we can study hadron 
production in p + p and A + A collisions at the LHC 
energies. At the highest energy of y/s = 14 TcV, mini- 
jet production involves partons in the very small x re- 
gion. The increase of the cut-off for independent mini-jet 
production required to fit the total cross section can be 
considered as indication of coherence or gluon saturation 
inside protons at very high energies. This will certainly 
affect the produced hadron pseudorapidity distribution 
at the LHC energies. Shown in Fig. @] is the energy de- 
pendence of the central pseudorapidity density of charged 



hadrons dN^/di] averaged over \rj\ < 1.0 as a function of 
the colliding energy in inelastic p + p collisions from both 
HIJING1.0 and HIJING 2.0 calculations as compared to 
experimental data, including that measured in p + p col- 
lisions at y/s = 0.9 TeV and 2.36 TeV from ALICE [U 
at the LHC. The central rapidity hadron density in HI- 
JING 2.0 continues to increase with the colliding energy 
at LHC more or less in a double logarithm form. Because 
of the larger values of cut-off for mini-jet production, the 
central rapidity hadron density in HIJING 2.0 is signifi- 
cantly lower than HIJING 1.0 at the LHC energies. 

Shown in Fig. [S] and [5] arc HIJING results for the pseu- 
dorapidity distributions of charged hadrons in p + p at 
= 0.9, 2.36 and 7 TcV in inelastic and non-single 
diffractive (NSD) events. The central rapidity density in 
NSD events is generally larger than in inelastic events 
and HIJING results agree with the experimental data 
reasonably well. Note that CMS 0, [47[ experiment uses 
a different definition of NSD events that gives slightly 
higher averaged multiplicity than that with UA5 defi- 
nition of NSD events in HIJING calculations. Another 
class of inelastic events in ALICE experiments (INEL>0) 
is defined by requiring at least one charged tracks within 
1 7/| < 1. HIJING results for these type of events (short- 
dashed lines) are similar to the inelastic events because 
the lack of double-diffractive events in the Mont Carlo 
model. Both single diffractive and non-diffractive events 
in HIJING have finite value of multiplicity in the central 
rapidity region. This disagrees with the ALICE data (4f| 
which are even larger than the central rapidity density 
in NSD events. Such discrepancy in trigger dependence 
of the averaged multiplicity can be fixed when double- 
diffractive events are included in future improvement of 
HIJING. However, this is not expected to have a signif- 
icant influence on particle production in heavy-ion colli- 
sions because the fraction of diffractive interactions be- 
comes increasingly smaller in multiple nucleon-nucleon 
collisions. 

Shown in Figs. [7] and [5] are the multiplicity distribu- 
tions for charged hadrons within different pseudorapidity 
range in inelastic and NSD p + p collisions at the LHC 
energies from HIJING calculations as compared to ex- 
perimental data from ALICE experiment [451 ]. HIJING 
results agree reasonably well with the data in low and 
intermediate multiplicities. However, they fall short of 
the experimental data at higher multiplicity tails in par- 
ticular at y/s = 7 TEV. These high multiplicity events 
are dominated by multiple jets and they are likely to 
have final state interaction with each other as indicated 
by the ridge structure in the two-hadron correlation in 
azimuthal angle and large rapidity gap observed in CMS 
experiment [481 ] . The final state interaction in these high 
multiplicity events could affect the multiplicity distribu- 
tion which is not included in current HIJING model. 

At high colliding energies, the transverse momentum of 
produced hadrons should be dominated by hadrons from 
jet fragmentation, especially at large transverse momen- 
tum, as shown in Fig. |H1 We compare the HIJING 2.0 
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FIG. 4: (color online) Central pseudorapidity density at r\ = 
for charged particles produced in p + p/p collisions as a 
function of y/s. The experimental data are from [§3, 52~5H . 
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FIG. 5: (color online) Pseudorapidity distribution dN c h/dr] 
in inelastic p + p collisions at LHC energies. The experiment 
data are from UA5[a3l, ALICEjH] and CMSEl, Sit Shown 
in short-dashed lines are the HIJING results for ALICE de- 
fined INEL>0 events in which at least one charged tracks are 
required within \rj\ < 1. 



results to the ALICE [49| experimental data on charged 
hadrons within \r)\ < 0.8 in inelastic p + p collisions at 
y/s = 0.9 TeV and CMS [H,[47| data on charged hadrons 
within |?7| < 2.4 in NSD p + p collisions at y/s = 0.9, 2.36 
and 7 TeV. The HIJING 2.0 results at y/s = 14 TeV 
for inelastic events in p + p collisions are also shown. 
Within the LHC energy range y/s = 0.9 - 14 TeV, the 
transverse momentum spectra at large pt region becomes 
apparently harder at higher colliding energies. Even 
though the cut-off for mini-jet production in HIJING 2.0 
is significantly larger than in HIJING 1.0, hadronization 
from partons as results of initial and final state radia- 
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FIG. 6: Pseudorapidity distribution dN c h/dr] in non-single 
diffractive (NSD) p + p collisions at LHC energies. The ex- 
periment data are from UA5[13, ALICEgl and CMS [H 53]. 
The definitions of NSD events in our calculations and the ex- 
perimental data from ALICE and UA5 follow that in UA5 
experiment as described in the text. CMS [HI, 53 uses a 
different definition. 
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FIG. 7: Multiplicity distributions for changed hadrons within 
\r)\ < 1 in inelastic p + p collisions at LHC energies from 
ALICE[45|] experiment compared to HIJING results. 



tion still fills in the hadron transverse momentum spec- 
tra in the intermediate transverse momentum. HIJING 
2.0 describes reasonably well the ALICE and CMS ex- 
perimental data at the LHC energies with the maximum 
discrepancy of 50% at small px ~ 2 GeV, as shown in 
the lower panels in Fig. [5] where the ratios of data over 
HIJING results are plotted. The transverse momentum 
spectra in this small px region can also be influenced by 
final-state interaction among jets. We also plot the ratios 
of charged hadron spectra from NLO pQCD calculation 
USUI over that from HIJING. HIJING reproduces the 
NLO pQCD results very well over large pr range and at 
different colliding energies. 

To describe hadron production in A + A collisions, one 
should include both the nuclear modification of the par- 
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FIG. 8: Multiplicity distributions for charged hadrons within 
different pseudorapidity intervals in NSD events of p + p colli- 
sions at = 2.36 TeV from ALICE[45| experiment as com- 
pared to HIJING results. 

ton distribution functions and jet quenching in final state 
interaction. Jet quenching in general suppresses high 
transverse momentum hadrons [26|. Since the yields of 
these high transverse momentum partons are relatively 
small as compared to the bulk part of the initial parton 
production, its effect on the total hadron multiplicity in 
A + A collisions will be small as illustrated by other mod- 
els studies (52rlo1| of hadron production in heavy-ion col- 
lisions at RHIC, most of which have also neglected the 
effect of final state interaction. With two parameters (po 
and <J S oft) fixed from fit to p+p collisions, the only uncer- 
tainty for hadron multiplicity density in A + A collisions 
is the nuclear modification of parton distribution func- 
tions, gluon distributions in particular, at small x. We 
assume parton distributions in nuclei, 

tf(x,Q 2 )=AR£( x ,Q*)fX( x ,Q2), (12) 

are given by that in nucleon and the modification fac- 
tor Ra(x,Q 2 ). This nuclear modification has been stud- 
ied with data from deeply inelastic scattering (DIS) and 
Drell-Yan lepton pair production experiments. However, 
most of parameterizations [56l loTj do not have gluon 
shadowing strong enough to give the central multiplic- 
ity density as measured in RHIC experiments within the 
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FIG. 9: (color online) Transverse momentum distributions of 
charged hadron in p + p collisions. The experimental data are 
from ALICE [49| (open square) for inelastic p + p collisions 
and CMS [46| (open and closed circle) for NSD p+p collisions 
at yfs = 0.9 and 2.36 TeV. The histograms are results from 
HIJING 2.0. 

two component model as implemented in HIJING model. 
We will use the HIJING new parameterization [29| which 
introduced a strong impact-parameter dependence of the 
parton shadowing in order to fit the centrality depen- 
dence of the central multiplicity density at RHIC. The 
shadowing factor for quarks and gluons are, 

Rf(x,b) = 1.0 + 1.19 log 1/f k (x 3 - 1.2.T 2 + 0.21x) 

- s q (b) (A^ a -l)°- 6 0.-3.5y/x) 

x exp(-x 2 /0.01), (13) 
Rg{x,b) = 1.0 + 1.19 log 1/f l4 (x 3 - 1.2x 2 + 0.21x) 

- Sg(b) (A^ 3 - 1) 06 (1 - 1.5X 035 ) 

x exp(-a; 2 /0.004), (14) 
respectively. The impact-parameter dependence of the 
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shadowing is implemented through the parameters, 

s a (b) = s a ^(l-b 2 /R A ), (15) 

where Ra = 1.12 A 1 / 3 is the nuclear size. The value 
s q = 0.1 is fixed by the experimental data on DIS off 
nuclear targets [29[ and s g = 0.17 — 0.22 is chosen to 
fit the RHIC experimental using HIJING 2.0 model as 
shown in Fig. 1101 The form of the impact parameter 
dependence is chosen to give rise to the centrality de- 
pendence of the pesudorapidity multiplicity density per 
participant pair and the range of values of s g is allowed 
by experimental errors at the RHIC energies. With this 
parameterization of parton shadowing, one can calculate 
the pesudorapidity multiplicity density per participant 
pair in Pb + Pb collisions at the LHC energies as a func- 
tion of the number of participants shown in Fig. 1101 As 
a comparison we also show results from HIJING1.0 as 
dashed lines in Fig. [TDJ The range of the gluon shadow- 
ing parameter s g allowed at RHIC leads to much bigger 
uncertainties in the multiplicity density at the LHC en- 
ergies. With such strong gluon shadowing constrained 
by the RHIC data, the effective mini-jet cross section at 
the LHC energies is still much larger than at RHIC and 
therefore the minijet component leads to stronger cen- 
trality (Np ar t) dependence of the central pesudorapidity 
density of produced hadrons. This strong centrality de- 
pendence is very different from other model predictions 
[58( such as the color glass condensate (CGC) or the 
limiting fragmentation model. Therefore, the first ex- 
perimental data on the hadron production at LHC will 
shed light on the parton/hadron production mechanism 
in high-energy hadron and nuclear collisions. 

To constrain parton shadowing in nucleus, one in prin- 
ciple can also study hadron production in d+A collisions. 
Shown in the lower panel of Fig. [TT] are the pseudorapid- 
ity distributions of charged hadrons in d + Au collisions 
at the RHIC energy ^/^200 GeV for different ccntralities 
and minimum-bias events as compared to the STAR data 
[65| . Within the experimental errors, both HIJING2.0 
(solid bands) and HIJING1.0 can describe the data rea- 
sonably well. We have followed the definition for the 
selection of centrality in the STAR measurements. The 
solid bands in the HIJING2.0 calculation as the varia- 
tion of the gluon shadowing parameter s g constrained by 
RHIC data of Au + Au collisions are too small to be dis- 
cerned by the d + Au data. Similar situation occurs in 
d + Pb collisions at the LHC energy as shown in the up- 
per panel of Fig. [TT] The change of the pesudorapidity 
distribution due to variation of the gluon nuclear shad- 
owing parameter is much smaller than in central Pb + Pb 
collisions shown in Fig. [UJ The difference between HI- 
JING2.0 and HIJING1.0 in d + Pb collisions becomes 
significant only at the highest energy s/s = 5.5 TcV, re- 
flecting the p + p results in Fig. 0] 




FIG. 10: (color online) Pseudorapidity density per partici- 
pant pair as a function of the average number of participants 
in heavy- ion collisions at different energies from HIJING2.0 
calculations. Dashed lines are HIJING1.0 results. The ex- 
perimental data are from Ref. [5!|. Note a vertical off-set is 
applied to each curve for clear illustration. 



A. Conclusions 

We have updated HIJING Monte Carlo model with 
modern parton distribution functions for the nuclcons 
and the new set of parameters within the two-component 
model for mini-jet production in high-energy nucleon- 
nucleon collisions. Because of the large gluon distribu- 
tion at small x in the GRV [28[ parameterization of the 
nucleon's PDFs used in HIJING 2.0, one has to intro- 
duce an energy-dependent transverse momentum cut-off 
Po for the mini-jet production and the soft parton inter- 
action cross section a so ft in order to describe the energy- 
dependence of the total, inelastic cross sections and the 
central rapidity hadron density in high-energy p + pip) 
collisions. The updated HIJING 2.0 model is shown to 
describe the existing experimental data on hadron pro- 
duction from ISR energy up to Fermilab Tevatron en- 
ergy The HIJING 2.0 results are also shown to be in 
good agreement with the recently published hadron spec- 
tra in p + p collisions at the LHC energies (-^=0.9, 2.36 
TeV), except for events with INEL^O trigger because of 
the lack of double-diffractive events in the model of soft 
interaction in HIJING. We also give the HIJING 2.0 pre- 
dictions ioT p + p collisions at yfs = 7 and 14 TeV. With 
a model parameterization for nuclear modification of the 
parton distribution functions, we also give HIJING 2.0 
prediction of hadron multiplicity in central Pb + Pb and 
minimum-bias events of d + Pb collisions at the LHC en- 
ergies y/s = 2.75 and 5.5 TeV/n. 

This is the first step of the upgrade to the HIJING 
model. Jet quenching description in current HIJING 
model is very schematic. The next stage of upgrade of HI- 
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FIG. 11: (color online) Pseudorapidity distributions of 
charged hadrons in d + Au collisions at the RHIC energy 
y^s = 200 GeV/n (lower panel) and d + Pb collisions (upper 
panel) at the LHC energies. The solids lines are HIJING2.0 
results with the bands corresponding to the variation of gluon 
nuclear shadowing parameter s g = 0.19 — 0.24. Data at the 
RHIC energy are from STAR experiment [65l ]. 



JING will be focused on jet quenching in dense medium 
incorporating the most recent development [60H64T ] in the 
theory of parton propagation and multiple interaction in 
dense medium. 
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